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Abstract 

Dynamic nonlinear (cubic) susceptibility in quantum d-dimensional Ising 
spin glass with short-range interactions is investigated on the basis of quan- 
tum droplet model and quantum-mechanical nonlinear response theory. Non- 
linear response depends on the tunneling rate for a droplet which regulates 
the strength of quantum fluctuations. It shows a strong dependence on the 
distribution of droplet free energies and on the droplet length scale average. 
Comparison with recent experiments on quantum spin glasses like disordered 
dipolar quantum Ising magnet is discussed . 
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The dynamics of glassy systems is an attractive and rapidly developing field of 
physics [1-5]. Spin glasses and quantum spin glasses are very interesting systems 
for theoretical and experimental investigation of dynamic phenomena [2-4]. Two 
major theoretical description of spin glass phenomena have evolved over the past 
twenty years:the mean field theory based on the Parisi replica symmetry-breaking 
approach and the droplet model based on renormalization group arguments [1-6]. 
The two pictures give very different physical interpretation of observable spin-glass 
phenomena. In this paper we investigate theoretically nonlinear cubic dynamic sus- 
ceptibility as function of frequency and temperature in the Ising spin glass in a 
transverse field in terms of quantum droplet model at very low temperatures (quan- 
tum regime). The quantum phase transition is governed by quantum fluctuations 
of the system which may tunnel from one local minimum of the free energy to an- 
other; new physical effects such as quantum channel of relaxation appear. There 
are few theoretical studies on the nonlinear static response in quantum spin glasses 
[7-10] and almost no studies on the dynamic nonlinear response [3]. In [3] dynamic 
nonlinear response of a quantum spin glass was found frequency independent and 
nonsingular in quantum critical regime in contrast to the behavior in usual spin 
glass. There are experimental data on the nonlinear dynamic response in classical 
[11-13] and quantum [14] spin glasses investigated by Fourier-transform technique. 
The third-order nonlinear susceptibility is negative and diverges at an ordinary spin 
glass transition temperature T g from both the upper and the lower sides. But when 
x' 3 is measured by a finite probing frequency the response falls out equilibrium be- 
fore the transition temperature and does not diverge at T g . Then shows a 
maximum at T ~ Tf(uj) where Tf{uj) is the freezing temperature which is the upper 
bound on T g and T g = Tf(uj — > 0). Such behavior was observed, for example, for 
classical Ising spin glass Feo.5Mr 5 Ti0 3 [13]. W. Wu et al. [14] measured nonlinear 
susceptibility ^{uj^T) in quantum spin glass (the diluted dipolar-coupled Ising spin 
glass LiHo .i67Yo.833F4 in the transverse field) tuning transverse field V from the 
T = classical to the T = quantum limit. At mK temperatures they found a 
clear dynamic signature of the spin glass to paramagnet transition whether domi- 
nated by thermal or quantum fluctuations. In [14] it was shown that the x'3 depends 
on frequency for uj > 10Hz. However it depends very weakly on uj for uj < 10Hz. 
There is a crossover between high uj (uj - dependent) and low uj (uj - independent) 
behaviors. Nonlinear susceptibility contains a diverging component which domi- 
nates at T = 98mK, but disappears by 25mK. The x'A 00 ) does n °t diverge but 
shows a maximum at Tf(uj). x'^) measured at a higher temperature and lower 
transverse field has a larger maximum than x'^) measured at a lower temperature 
and larger transverse field. The analysis of these experimental data seems not clear 
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[15] because frequencies used in the experiments [14] are not sufficiently low such as 
to determine the equilibrium behavior of system. Contrary to the theoretical expec- 
tations, quantum transitions may be qualitatively different from thermally driven 
transitions in real spin glasses. Recently the linear dynamic susceptibility (in-phase 
and out of phase components) at T = was investigated theoretically for the Ising 
spin glass in a transverse field in terms of quantum droplet model by M. J. Thill 
and D. A. Huse [6]. A basic assumption of the droplet picture is that the spin glass 
dynamics is governed by large-scale excitations whose relaxation time increases with 
length scale. In previous papers [16] we have calculated the real and imaginary parts 
of linear dynamic susceptibility in the same model, as in [6], for very low nonzero 
temperatures using the general linear response theory of magnetic dispersion and 
absorption phenomena for quantum systems by Kubo and Tomita [17]. We note 
that in [6] the real part of the cubic nonlinear ac susceptibility was defined as the 
in-phase 3u magnetization response M(3uj) to a small time-dependent applied field 
hcos(ut) 

where V is the sample volume. The authors of [6] gave some expression for x'3 they 
only expect at zero temperature. 

The full nonlinear response theory despite its generality and importance is of 
limited practical value because it is mathematically difficult. It is necessary to make 
approximations, as the well-known perturbation expansion of the time-evolution 
operator, using Hamiltonian with some small parameter. Nonlinear response theory 
was developed and described, for example, in [17-22]. Here we summarize briefly 
the theory of higher-order dynamic response. It is based on the Hamiltonian 

Ht = Ho + Hi = H - AjFjd), t > t (2) 

where Ho is nonperturbation Hamiltonian of system, Hi is perturbation Hamilto- 
nian which describes the interaction between the Heisenberg operators Aj (material 
operators) and the external perturbation Fj. At time t > t one is interested in the 
expectation value of the Heisenberg operator Bi which is given by 

(Bi(t)) = Tr[poMt)] = Tr[p(t)Bi\ (3) 

where the density matrix p(t) = U(t,to)p(to)U^(t,to). The time-evolution operator 
U satisfies the Schrodinger equation ih^r = HU. It is difficult to find an expression 
for U in closed form. It was used that Hi(t) is in some sense small. We define the 
total response of the system at time t to the external force Fj as the difference 

ABi(t) = (Bi(t)) - (Bi) (4) 
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where the subscript zero on expectation values refers to the equilibrium expectations. 
One can then understand the behavior of the system in terms of the dynamical 
response. Using aforementioned expressions [2-3] the dynamical response can be 
written through the third order in the perturbation 7Yi in the following form 

rt 



(B i {t))-(B i ) Q = &B i {t)~ I dt' lfij (t-t')F j (t)+ 

J to 

f dt x t dt 2 ip ijk (t -t 2 ,h- t 2 )F k (t 1 )F j (t 2 )+ 

J to J to 

f dh f 1 dt 2 f 2 dh (fiijUt - h, h - h, h - ^FMF^Fjfo) + ... (5) 

J to J to J to 

where <fij, iptjk, tfiijki are the first-, second- and third order response functions, 

<p ij (t-t') = ±([A j ,B i (t-t')]) , (6) 



<p ij k(t-t 2 ,t 1 -t 2 ) = ^{[Aj^Akfa-hlBiit-tzWo, (7) 



<PiMt - h, t 2 - 1 3 , h - t 2 ) = jr^{[Aj, [A k (h - t 2 ), [Mh - 1 3 ), Bi(t - t 3 )]}}) . (8) 

Here we have employed the summation convention over repeated indices and cyclic 
invariance of the trace. The expressions (6-8) can be written in a more revealing 
form if we set to = —oo and change integration variables including an adiabatic 
switching factor if necessary. We take (instead of t, t±, t 2 , t 3 ) T\ — t — t\, r 2 — t\ — t 2 
time differences further. Ordinary linear response theory utilized only </ 7 ij( T 'i) an d 
it is the simplest approximation to the full theory of linear dynamic response [23]. 
Using new variables we may write the expressions for response functions in the form 

^—^(^(O)]), (9) 

^ k = j^ 2 (U{n + T 2 lA{T 2 )],B\), (10) 

<P% = -^([[[Mn + T2 + r 3 ), A(r 2 + r 3 )], A(r 3 )], B\) (11) 

where the bracket (...) denotes an expectation value with respect to the equilibrium 
ensemble. The useful interpretation is generated from eq. (5) in the case that 
to = — oo if we consider the external force F is constant and vanishes for t > 0. For 
t — the system is in partial equilibrium and starts to relax to equilibrium. It is 
convenient to write nonlinear response for this case (initial value case [21]) as 



4 



(B(t)) - (B)q = R {1 \t)F + \R {2 \t)FF + \R$\t)FFF + ... (12) 
where R a (t) are the relaxation functions, 

Rf»(t) = fdr^fin), (13) 

R (2 \t) = j o dnj o dr 2 ^%(r 2 ,T 1 + r 2 ), (14) 

roc /■oo /■oo , , 

R®(t)=J Q dnj^ dr 2 j o rfr 3 4 3 Mr3,r 2 + r3,r 1 + r 2 + r3). (15) 



In this form response may describe relaxation of the system. If response function 
(f(t)^ vanishes as t — > oo, then ip(t)^ = —^R^(t), so R^(t) contains much more 
information than the response function. 

Let the ac magnetic field h w = hcos(cut) be applied to a magnetic system. The 
magnetization nonlinear response M{u,t) = Sfc=i{@'fe cos(kut) + 0^ sin(ku;t)} to 
harmonic magnetic field h contains only odd harmonic; X ~ x'iA © 3 ~ Xs^ 3 , 
etc [11]. In expression for M(u,t) the magnitudes Q' k and 0j£ are real and imagi- 
nary parts of harmonic amplitudes respectively. For the general theory of nonlinear 
processes one can evaluate [11] 

h 3 

©1 =XiM)H [x' 3 {uj,0,uj)]—+ [Ax' 5 {uj,0,uj,0,uj) + 2x' 5 {uj,0,u;,2uj,uj) + 



2x' 5 (w, 2 ^> ^> °> + Xe>( w » 2u;, w > 2w, w) + Xs(^, 2w, 3w, 2u, + (16) 

@3 = X 3 (3^,2w,w)— + [X5(3w,2w,3w,2^,w) + X5(3^,4^,3w,2w,w) 

2Xs( 3 ^> 2w, 0, w) + Xs(3^, 2w, w)] — + ... , (17) 

©; = X5 (5^,4^,3^,2^)— + ... (18) 



The measurement of all the harmonic amplitudes 0^ gives a measurement of the 
susceptibilities Xa in two limits: a) if X\h ^> X^h 3 ^> Xs^ 5 the back reaction is 
negligible and each harmonic measures the susceptibility of the same order; b) in 
the static (u — > 0) limit the solution of the linear system [15] fully accounts for the 
back reaction. In the absence of ac magnetic field, the back reaction can be made 
small so that the dynamic susceptibilities can be obtained from eq. (15). In a more 
compact notation we may write 

M(u, t) ~ [M + M U) + M 3uJ + ...] + (complex conjugate) (19) 

where M is the equilibrium magnetization in zero field; M w is the cu-magnetization 
response; M 3u) is the 3a;-magnetization response and so on. 
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The expressions (9-11) may be considered as solution of the corresponding quan- 
tum equations considered above. The external ac field we assume classical value. 
This field interacts with quantum system and system behavior is determined by 
quantum laws. We shall focus on the real part of the third-order nonlinear dynamic 
susceptibility x' 3 (3cj, 2cu, and denote it as ^(u). In this paper we are interested 
in the response when the ac magnetic field is applied in z-direction; xi — Xzz and so 
on. In formulas (2-11) both Ai and Bi are the magnetic dipole moment operators. 
Considering the initial value case [21] we suppose, like in [6], that the system is in 
equilibrium with a small time-independent field h — h z for t < and that the exter- 
nal field is turned off at t = 0, then for t > the induced magnetization of the sample 
in z-direction to first order of perturbation theory is given by M(t)-M « RM(t)h, 
where a relaxation function 

RW(t) = / ^){ T )dT (20) 
Jo 

and the first order response function is [20] 

^{r)^-^m{rlM 3 ]) . (21) 
The higher-order response functions are given by 

4 2 ^i,r 2 ) « ^([[M,(r 1 + r 2 ),M J (r 2 )],M fc ]) ; (22) 
4 3 L(ri,r 2 ,r 3 ) « -JLfl^fa + r 2 + r s ), Mj(t 2 + r 3 )] , M fc (r 3 )] , M l \). (23) 

The linear and the nonlinear dynamic susceptibilities (admittances in the spectral 
representation [22]) may be found through response functions (21-23). In order 
to find the complete expression for susceptibility, we should use its symmetry and 
causality properties [22]. 

The nonlinear susceptibilities may be chosen so that these susceptibilities were 
symmetrical relative to simultaneous permutation of tensor indices and correspond- 
ing to them arguments, for example, the second rank tensors are Xijk(wi,u)2) = 
Xikj(w2,Wi), and the fourth rank tensors are 

Xijkifai, ^3) = Xikji(u 2 , wi, u 3 ) = Xijiki^u ws, cj 2 ) = ... (24) 

according to causality property 

Xij = for ti < max (r 2 , r 3 , . . .). (25) 

Linear and nonlinear dynamic susceptibilities are given by 

X«H = / o °° dr^\r)e^ (26) 
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' P OO P OO 

Xvkfau*) = yj o dnJ Q rfr 2 {4 2 i(r 1 ,r 2 )e^ + -)- + — } (27) 

Xijfci(wi,w2,w3) = y o <foi y o rfr 2 y o dT 3 <p\ j > kl {T 1 ,T 2 ,T 3 ) x 

x exp [i(cji + cj 2 + ^3)ti + 2^ 2 r 2 + icj 3 r 3 ] (28) 

In eq. (28) P 3 means the summation over all permutations of the subscripts (ouij), 
(u; 2 /c) and (cu^l) [22]; response functions ip are given by expressions (21-23). 

In particular, the second harmonics generation process is characterized by tensor 
Xijk(u,u). If uj\ = u)2 = UJ3 = cu, the triple frequency 3u is formed; the frequency 
tripling is described by tensor Xijki(w,u),w) (we note it as x 3 (o;)). 

The droplet model of classical Ising spin glass was considered by D. S. Fisher 
and D. A. Huse [4]. The features of this model are described also, for example, in 
[23]. In the droplet model there are only two pure thermodynamical states related 
to each other by a global spin flip. In magnetic field there is no phase transition. 
A droplet is an excited cluster in an ordered state where all the spins are inverted. 
The natural scaling ansatz for droplet free energy e L (which are considered to be 
independent random variables) is ~ L e ,L > ((T); ( is the correlation length, 
L is the length scale of droplet and 9 is the zero temperature thermal exponent, 
9 < (d — l)/2. One droplet consists of order L d spins. Below the lower critical 
dimension d t , 9 < 0; above d t one has 9 > 0. 

Recently M. J. Thill and D. A. Huse [6] have shown that the rf-dimensional 
quantum ising spin glass in a transverse field with Hamiltonian 

w = -XX-s?s;-r£5? (29) 

(where Si are the Pauli matrices, T is the strength of the transverse field and the 
nearest neighbor interactions 1^ are independent random variables of mean zero) can 
be represented as the Hamiltonian of the independent quantum two-level systems 
(low energy droplets) of the form 

"^EE^A + r^J (so) 

L D L 

where S z Dh and Sf) are the Pauli matrices representing the two states of the droplet; 
the sum is over all droplets D L at length scale L and over all length scales L, and 

~ r°° dL 

?~1t < 31 > 

with a short-distance cutoff L . The value T L regulates the strength of quantum 
fluctuations (T L — > corresponds to the classical limit). 

r £ = r e- CTLd (32) 
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is the tunneling rate for a droplet of linear size L, Tq is the microscopic tunneling 
rate; a is defined from the equation 2K = aL d where 2K is the surface free energy 
of an interface between the two droplet states, so a is a reduced surface tension for 
this interface; a is approximately the same for all droplet. We will assume is the 
same for all droplets of scale L. The droplet excitations have a broad distribution 
of their free energies at scale L for large L in a scaling form [4,6] 

p ^ = 4w v fe) ' L ^ 00 ■ (33) 

It is assumed that P L (x -> 0) > 0, P L (0) - P L (x) ~ at x -> 0. 7(T) is a gen- 
eralized temperature dependent stiffness modulus which is of order of characteristic 
1 

exchange X = (Zfj) 2 at T = and vanishes for T >T g . 

There is a crossover length scale, L*(T), defined by condition r^*^) = k B T or 
L*(T) = (ilog^) 1/d . For droplets with L < L*(T) and Y L > k B T the energy 



e l + r| i s always more than k B T and thermal fluctuations are insignificant at 
temperature T. Droplets with L ^> L*(T) have I\ A;^T and behave classically. 
The large droplets (e^ < k B T, < k B T) are thermally active. At low T only a 
small fraction of droplets is thermally active, but many low-T static properties are 
dominated by these droplets at the crossover length L*(T). 

The total magnet moment of a droplet will scale as qL% where q a random 
number with mean zero and q 2 ~ q B A, qEA = (Sf) 2 is the Edwards-Anderson order 
parameter [4]. The total magnetization M of the sample will be 

^ = E^E(^)^ f (34) 



L L '' D L 



and 

M 



m = y = Y: (S' DL )qL* (35) 

L 



where (S z DL )qL~ means the average over the droplets energies e L . The static sus- 
ceptibilities are defined in terms of the expansion of magnetization into Taylor series 
approximately as 

M = xih-xsh 3 + ... . (36) 

For enough small external field it is possible to be restricted by few terms of this 
expansion. The susceptibilities are got as derivatives with respect to h at h — : 



X* 1 dh 



h=0 



X3 - QhS- 



h=0 



M. J. Thill and D. A. Huse have calculated static linear and nonlinear suscep- 
tibilities for droplet system described by the Hamiltonian (29). The static linear 
susceptibility diverges at T = below the lower critical dimension di. The static 
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nonlinear susceptibility diverges in all dimensions d. The static linear susceptibility 
appears to start away from the nonzero constant T = value decreasingly versus T 
to lowest order [6]. 

Now using the quantum droplet model of the short-range Ising spin glass in 
a transverse field and quantum- mechanical case (3Tl ^> 1 (quantum regime) we 
calculate the third order nonlinear dynamic response at very low finite temperatures. 

When we consider the droplets at finite temperatures they may have two char- 
acteristic rates, (a) the Rabi frequency (is of order Ti) and (b) the rate of classical 
activation over energy barrier B for annihilation and creation of the droplet excita- 
tions [6] 

t ~ r e^ , (37) 

where B ~ AL^, < ip < d — 1, ^ is some exponent [4], A is a barrier energy 
at T < T 9 , A ~ X; r is a microscopic time. There is a complicated dynamical 
classical-to-quantum crossover depending on temperature, frequency of ac external 
field and length scale L. According to [6] the crossover dynamic length is determined 
from the condition T^ 1 = T, i.e. 

L* dyn {T) - {^k B T)^ d (38) 

The system behaves presumably classically or quantum mechanically when the dom- 
inant length scale L is above or below L* dyn for frequency u. 

Now we consider dynamic third-order susceptibility x'z^iT) at finite very low 
temperatures (quantum regime) when f3T L ^> 1. We define nonlinear third order 
dynamic susceptibility x'zd^iT) by expression (23) and (28) 

X 3 (^> T ) = jj^ J o dl ~i J o dr 2 J Q dr z exp [i (Sun + ^t 2 + ujt 3 )] x 

([[[Miin + r 2 + r 3 ), Mj(T 2 + r 3 )] , M k (r 3 )) , Mi]). (39) 

The contribution of a single droplet to the real part of dynamic third-order suscep- 
tibility up to some factor ~ q 2 EA L 2d is proportional to 

L 



(ELo A k (u, T L )ef) tanh [\^e 2 L + Tj ^ 
~ ^ (el + r|)f (el + Tl - 9u 2 ) (el + T 2 L - u 2 f (e| + T| - £ ) " ' ^ 



where A = u 2 ^ +4^1-5^1; A l = 2r^-12, 5u 2 T 8 L -49, 75cu 4 r|-23, 5o; 6 r|- 
2, 25u 8 T 2 L ; A 3 = 12T 6 L - 45, 5u 2 T 4 L + 41, 75u 2 T 4 L ; A A = 8T 4 L - 15, hu 2 Y 2 L ; A 5 = 2T 2 L . 

We have to average X3D L over droplet energies using the distribution of droplet 
free energies (33) and changing variables from to x = fie^. 
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After averaging over droplet energies for cases Tl > 3lo,Tl ~ 3uj,Tl < 3lo we 
receive that the real part of the nonlinear susceptibility is dominated by droplets of 
length scale 

1 



J dom(3ui) 



(41) 



which is determined by condition T L ~ 3uo. Then for T L > 3uj the result of two 
averages is given by the following expression 



X 3 ~ 
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|tt sec 
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^-fe K0-fe))" 

d 



-G 



a. 



(<f>-k) 



+ 



ad 



0-5-2fc 



(42) 



fe=0,l 



where G[a, z] is incomplete gamma-function. The coefficients in this expression 
depend on only and are given by 



a = 



d-6(l+<t>) . 



4 o — — i — ^ 4 



B = G 
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fc=l 
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for T L ~ 

C = 
C 2 = G 
C 3 = -G 
C4 — — G 

This expression does not diverge if 2 < < 3 and —1 + 4 < 0. This expression 
has singularity at u ~ When the frequency increases the values of X3 are growing 
to infinity while u — > X3 maintains this property when a; is more than The 
temperature dependence of X3 in this case has no extremes, we observe monotonous 
decrease of values of x'3 with temperature. 

Nonlinear susceptibility given by expression (42) does not diverge if 2 < < 3 
and -1 + I < 0. 

We observe that nonlinear susceptibility has strong dependence on distribution 
function Pl(^l), i- e. on 0, on droplet microscopic tunneling rate To and other 
parameters. One can see that the real part ^ 3 {u,T) varies approximately logarith- 
mically with frequency. This signalizes broad distribution of relaxation times of the 
system. 
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Let us take for numerical calculation the following numbers: 9 = 1,Tq = 
10 10 s _1 ,ci = 3, = 2.5,7 = 10~ 15 erg,a = 10~ 15 ,g = 0.5. The frequency depen- 
dence of x^i^^T) is shown in Fig. 1. We give log 10 /-dependence at log 10 / from 
to 11 at fixed several temperatures: T\ = 0.001, T2 = 0.005, T3 = 0.01, T4 = 0.05. 

The frequency interval covers some decades of frequencies. Our numerical calcu- 
lations show the crossover between Iow-cj and high-u; behaviors. In low-u; region the 
nonlinear response is found nonsingular and slowly decreasing. When frequency in- 
creases the curve falls down more quickly, the nonlinear response diverges at cu ~ 
then the curve rises to some value. In low-u; region we have a qualitative agreement 
with experimental data for disordered dipolar magnet LiHo :r Y 1 _ ;I .F4. At different low 
fixed temperatures the behavior of x^i^^T) is the same and the values of rf 3 (iv,T) 
are approximately the same. Therefore we give only one curve for all fixed temper- 
atures. 

In Fig. 2 we give the temperature dependence of x 3 (u;,T) at temperatures from 
to 10~ 2 K at fixed several frequencies: fx = 10 7 Hz, f 2 = 2.5 x 10 7 Hz, f 3 = 
5 x WHzJi = 7.5 x 10 7 Hz, f 5 = 10 s Hz of ac field (/ = §). The behavior 
of ^3(0;, T) indicates the following glassy-like features. The curves of the temper- 
ature dependence of x'^if^iT) have maxima depending on fixed frequency. The 
temperature of ^-maximum Tf{uj) depends on frequency. The nonlinear suscepti- 
bility magnitudes at different fixed frequencies are remarkably distinguishable. The 
temperatures of maximum values are different. When the frequency increases the 
temperature of x 3 -maximum shifts towards high temperatures. The similar curve of 
temperature variation of x'3 was observed in spin glasses at more high temperatures 
[2,13]. If we consider only T-dependent part of x'3 we see that the ^-maxima are 
sharp (Fig. 3). 

The cubic dynamic susceptibility x 3 (3o;) is analitically and numerically calcu- 
lated in quantum spin glass in terms of quantum droplet model on the basis of 
general dynamic nonlinear quantum-mechanical response theory. We have carefully 
analyzed the susceptibility temperature-frequency behavior to study the properties 
of the low temperature magnetic state and to determine whether or not a conven- 
tional spin glass state exists below Tf. Comparing with the case of a true spin glass 
transition we see that our data indicate that the magnetic state below Tf does not 
correspond to a conventional spin glass state below Tf. We find a glassy type slow 
dynamics. Similar frequency dependence was observed by W. Wu et al.[14]. 

Our calculations at T = coincide with T = result of M. J. Thill and D. A. Huse 
[6]. For finite temperatures we find some features which have been recently observed 
[13]. We suppose that at some very low temperature Tf (temperature of maximum 
of Xsi^j T)) there is a phase transition. If 6 > and d = 3 we suppose a true phase 
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transition at very low temperature Tf ~ 1CT 4 -j- 8.5 x 10~ 4 K for / = 10 7 -j- 10 8 Hz 
respectively (Fig. 3). 

Besides frequency and temperature dependence the shape of Xs(3w) depends 
crucially on the probability distribution of droplet free energies, on the tunneling 
rate for a droplet of linear size L, on the material parameters. In consequence of 
this dependence there is divergence (or convergence) of x 3 (3o;). We need to take 
into account (in future paper) the dipole-dipole interaction between droplets and 
also droplet-lattice interaction. 

Applying our results to the reported experimental data on the nonlinear dy- 
namic susceptibility of LiHo a; Y 1 _ 2 .F4 we observe that a fairly good agreement may 
be achieved. 
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Figure captions 

Fig.l - The frequency dependence of the real part of the nonlinear dynamical 
susceptibility at fixed temperature. 

Fig. 2 - The temperature dependence of the real part of the nonlinear dynamical 
susceptibility at various frequencies. 

Fig. 3 - The temperature dependence of the T-dependent part of the nonlinear 
dynamical susceptibility at various frequencies f. 



15 



This figure "Figurel.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/020350 1 v 1 



This figure "Figure2.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/020350 1 v 1 



This figure "Figure3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/020350 1 v 1 



